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Abstract. The most detailed existing proposal for the structure of 
spacetime singularities originates in the work of Belinskii, Khalatnikov 
and Lifshitz. We show rigorously the correctness of this proposal in 
the case of analytic solutions of the Einstein equations coupled to a 
scalar field or stiff fluid. More specifically, we prove the existence of 
a family of spacetimes depending on the same number of free func- 
tions as the general solution which have the asymptotics suggested by 
the Belinskii-Khalatnikov-Lifshitz proposal near their singularities. In 
these spacetimes a neighbourhood of the singularity can be covered by 
a Gaussian coordinate system in which the singularity is simultaneous 
and the evolution at different spatial points decouples. 
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1. Introduction 

The singularity theorems of Penrose and Hawking are among the best 
known theoretical results in general relativity. They guarantee the existence 
of spacetime singularities under rather general circumstances but say little 
about the structure of the singularities they predict. In the literature there 
are heuristic approaches to describing the structure of singularities, notably 



that of Belinskii, Khalatnikov and Lifshitz (BKL), described in [14], J4J], || 
and their references. The BKL work indicates that generic singularities are 
oscillatory and therefore, in a certain sense, complicated. This complexity 
may explain why it has not been possible to determine the structure of the 
singularities by rigorous mathematical arguments. 
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According to the BKL analysis, the presence of oscillatory behaviour in 
solutions of the Einstein equations coupled to some matter fields is to a 
great extent independent of the details of the matter content. There are, 
however, exceptions. It was pointed out by Belinskii and Khalatnikov^] 
that a massless scalar field can change the situation dramatically, producing 
singularities without oscillations. A massless scalar field is closely related 
to a stiff fluid, i.e. a perfect fluid with pressure equal to energy density, as 
will be explained in more detail below. Barrow Q exploited the singularity 
structure of solutions of the Einstein equations coupled to a stiff fluid for 
a description of the early universe he called 'quiescent cosmology'. We will 
refer to singularities where oscillatory behaviour is absent due to the matter 
content of spacetime as quiescent singularities. 

Recently the Einstein equations coupled to a scalar field have once again 
been a source of interest, this time in the context of string cosmology. A 
formal low energy limit of string theory gives rise to the Einstein equations 
coupled to various matter fields. Under simplifying assumptions the col- 
lection of matter fields can be reduced to a single scalar field, the dilaton. 
The field equations are then equivalent to the standard Einstein-scalar field 
equations. (Note, however, that the metric occurring in this formulation of 
the equations is not the physical metric.) The structure of the singularity 
in these models plays a role in the so-called pre-big bang scenario. For more 
information on these matters the reader is referred to the work of Buonanno, 
D amour and VenezianojT]]. 

In view of the above facts, the Einstein-scalar field equations and, more 
generally, the Einstein-stiff fluid equations represent an opportunity to prove 
something about the structure of spacetime singularities in a context simpler 
than that encountered in the case of the vacuum Einstein equations or the 
Einstein equations coupled to a perfect fluid with a softer equation of state. 
In this paper we take this opportunity and prove the existence of a family 
of solutions of the Einstein-scalar field equations whose singularities can be 
described in detail and are quiescent. These spacetimes are very general in 
the sense that no symmetry is assumed and they depend on as many free 
functions as the general solution of the Einstein-scalar field equations. They 
have an initial singularity near which they can be approximated by solutions 
of a simpler system of differential equations, the velocity dominated system. 
Like the full system, it consists of constraints and evolution equations. The 
evolution equations contain no spatial derivatives and are thus a system of 
ordinary differential equations. This is an expression of the idea of BKL 
that the evolution at different spatial points decouples near the singularity. 

The structure of the paper is as follows. In the second section we recall 
the Einstein-scalar field and Einstein-stiff fluid equations and define the 
corresponding velocity dominated systems. This allows the main theorems 
to be stated. They assert the existence of a unique solution of the Einstein- 
scalar field equations or Einstein-stiff fluid equations asymptotic to a given 
solution of the velocity dominated system. The proofs of the existence and 
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uniqueness theorems are described in the third section. In the fourth section 
the main analytical tool used in these proofs, the theory of Fuchsian systems, 
is presented. The algebraic machinery needed for the application of the 
Fuchsian theory is set up in the fifth section. This provides the basis for 
the estimates of spatial curvature and other important quantities in the 
section which follows. The seventh section treats relevant aspects of the 
constraints. The paper concludes with a discussion of what can be learned 
from the results of the paper and what generalizations are desirable. 

Throughout the paper the scalar field and stiff fluid cases are treated in 
parallel. These are independent except in section [/] where the propagation of 
constraints for the scalar field is deduced from the corresponding statement 
for the stiff fluid. Hence concentrating on the scalar field case on a first 
reading would give a good idea of the main features of the proofs. 

2. The main results 

Let ^g a /3 be a Lorentz metric on a four-dimensional manifold M which is 
diffeomorphic to (0, T) x S for a three-dimensional manifold S. Let a point 
of M be denoted by (t, x), where t € (0, T) and x £ S. It will be important 
in the following to express the geometrical quantities of interest in terms 
of a local frame {e a } on S. Let {9 a } denote the coframe dual to {e a }. 
Throughout the paper lower case Latin indices refer to components in this 
frame, except where other conventions are introduced explicitly. Suppose 
that the metric takes the form: 

(2.i) - dt * + g ah (t)e a ®e b 

where g a bif) denotes the one-parameter family of Riemannian metrics on 
S defined by the metrics induced on the hypersurfaces t =constant by the 
metric ( 4 ) g a p . A function t such that the metric takes the form ( |2.lD is called 
a Gaussian time coordinate. In this case, the second fundamental form of a 
hypersurface t =constant is given by k ab = — \dtg a b- 

2.1. The Einstein-matter equations. The Einstein field equations cou- 
pled to matter can be written in the following equivalent 3 + 1 form. The 
constraints are: 

(2.2a) R - k ab k ab + (tr/c) 2 = 16np 

(2.2b) V a k ab - e b (tvk) = 8irj b 

The evolution equations are: 
(2.3a) d t g a b = -2k ab 

(2.3b) d t k\ = R a b + (trk)k a b - 8ir(S a b - ^trS) - 4vr^ a 6 

Here R is the scalar curvature of g ab and R ab its Ricci tensor. The quantities 
p, j a and S ab are projections of the energy-momentum tensor. Their explicit 
forms in the cases of interest in this paper will be given below. 
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The energy-momentum tensor of a scalar field is given by 

(2.4) T ap = V a 4>Vp4> - ^(V 7 </>V^) (4) <fo/3 

The Einstein equations can be written in the equivalent form ^ R a p = 
8irV a 4>V p4>, where ^ ' R a p is the Ricci tensor of ^g a p- We have p = Tqq, 
j a = —To a and S a b = T ab so that in the case of a scalar field it follows from 
(HD that: 

(2.5a) p=^[(d t ct ) ) 2 + g ab e a ^)e b ^)} 

(2.5b) j b = -d t 4>e b {4>) 

(2.5c) S ab = e a (0)e fe (0) + ^[{d^f - g cd e c (0)e d (<j))}9ab 

Note that it follows from the Einstein-scalar field equations as a consequence 
of the Bianchi identity that <fi satisfies the wave equation ^ g a/3 V a V p<j) = 0. 
This has the 3 + 1 form: 

(2.6) -df(p + (\xk)dt<j> + Acf) = 

The constraints and evolution equations are together equivalent to the full 
Einstein-scalar field equations. In the following we will work with the 3 + 1 
formulation of the equations rather than the four-dimensional formulation. 

A stiff fluid is a perfect fluid with pressure equal to energy density. As 
we will see, it is closely related to the scalar field. The energy-momentum 
tensor of a stiff fluid is 

(2.7) T a/3 = p(2u a U[3 + ^g af 3) 

where p is the energy density of the fluid in a comoving frame and u a is the 
four-velocity. The Euler equations are obtained by substituting this expres- 
sion into the equation V a T a ^ = 0. The relation between the scalar field and 
the stiff fluid is as follows. Given a solution of the Einstein equations coupled 
to a scalar field, where the gradient of the scalar field <f> is everywhere time- 
like, define fi = -(l/2)V a 0V a and u a = ±(-V /3 (/)V /3 0)~ 1 / 2 V a (/>. Here the 
sign is chosen so that u a is future pointing, and so can be interpreted as the 
four-velocity of a fluid. Then the energy-momentum tensor defined by (|2.7| ) 
is equal to the energy-momentum tensor of the scalar field. Since the latter 
is divergence- free we see that the fluid variables just defined together with 
the original metric define a solution of the Einstein equations coupled to a 
stiff fluid. In the spacetimes of interest in the following, the gradient of <p is 
always timelike near the singularity, so that the condition on the gradient is 
not a restriction in that situation. 
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The matter terms needed for the Einstein equations are given in the stiff 
fluid case by 

(2.8a) p = ^(l + 2|u| 2 ) 

(2.8b) j b = 2 M (1 + \u\ 2 ) l l 2 u b 

(2.8c) S ab = fi(2u a u b + g ab ) 

Here \u\ 2 = g ab u a u b . The Euler equations can be written in the following 
3+1 form: 

dt[i — 2(tik)n = — 2\u\ 2 dtfJ, ~ 4/ro a "dtu a — 2fik ab u a u b 
+2(trk)n\u\ 2 - 2e a (n)(l + \u\ 2 ) l ' 2 u a - + \u\ 2 )~ 1/2 u a V a u b u b 

(2.9a) -2/i(l + |n| 2 ) 1/2 V a n a 

d t u a + {txk)u a = -^[dtn - 2(tik)fj]u a + (1 + j-ul 2 ) -1 ^^^^ 

(2.9b) -(1 + \u\ 2 )- l l 2 [{u a u b + (l/2)^) / u~ 1 e 6 (^) + (V c u c u a + u c V c u a )] 



2.2. The velocity dominated system. We will prove the existence of a 
large class of solutions of the Einstein-scalar field equations and Einstein-stiff 
fluid equations whose singularities we can describe in great detail. These 



singularities are of the type known as velocity dominated. (Cf. || and [11].) 
This means that near the singularity the solution can be approximated by 
a solution of a simpler system, the velocity dominated system. Like the 
3+1 version of the full equations it consists of constraints and evolution 
equations. Solutions of the velocity dominated system will always be written 
with a left superscript zero and the convention is adopted that the indices of 
all quantities with this superscript are moved with the velocity dominated 
metric °g c d- The velocity dominated equations will now be written out 
explicitly. The constraints are: 

(2.10a) -°k ab °k ab + (tr°A;) 2 = I6ir°p 

(2.10b) V a {°k ab ) - e b (ti°k) = 8ir°j b 

The evolution equations are: 

(2.11a) d t °g ab = -2°k ab 

(2.11b) d t °k a b = (tr°k)°k a b - 8Tr(°S a b - ^6%tr°S) - 4iv p6 a b 



In these equations the matter terms are not identical to those in the full 
equations but have been obtained from those by discarding certain terms. 
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In the case of a scalar field the (truncated) energy-momentum tensor 
components are given by 

(2.12a) ° p =I(d t » 2 
(2.12b) °j b = -d t V b 0V) 

(2.12c) °S ah = ^(d t <P) 2 (°g ah ) 

The scalar field satisfies the equation 

-d t 2 (» + (tr fc)d t V = 

It is important to note that the velocity dominated evolution equations are 
ordinary differential equations. However the constraints still include partial 
differential equations. 

In the stiff fluid case the (truncated) energy-momentum tensor compo- 
nents are 

(2.13a) °p = V 

(2.13b) °j b = 2° Aft 

(2.13c) °S ab = W 

The velocity dominated Euler equations are 

(2.14a) dt°n - 2(tr°fc) V = 

(2.14b) d t °u a + (tv°k)°u a = -(l/2)V _1 e a (V) 

In contrast to the scalar field case, this is not a system of ordinary differential 
equations. However it has a hierarchical ODE structure in the sense that if 
the ODE for V is solved and the result substituted into the other equations 
an ODE system for the °u a results. 

Substituting the expressions for the truncated energy-momentum tensor 
into the velocity dominated system shows that the matter terms in the 
velocity dominated evolution equation for °k ab cancel both for the scalar 
field and the stiff fluid, leaving 

d t °k a b = (ti°k)°k a b 

Taking the trace of this equation gives dt(ti°k) = (tr°/c) 2 . This has the 
general solution tr°/c = (C — If we wish an initial singularity, as 

signalled by the blow-up of tr°/c, to occur at t = then tr°/c = — t -1 . Going 
back to the equation for °k a b we see that t(°k a b ) is independent of time. 
Thus all components of the mixed form of the second fundamental form are 
proportional to t -1 . 

At any given spatial point we can simultaneously diagonalize °g ab and 
°k ab by a suitable choice of frame. The matrix of components of the metric 
in this frame is diagonal with the diagonal elements being proportional to 
powers of t. This form of the metric is that originally used by BKL. Its dis- 
advantage is that in general this frame cannot be chosen to depend smoothly 
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on the spatial point. (There are difficulties when there are changes in the 
multiplicity of the eigenvalues of °k a b .) This is one reason why a different 
formulation is used in this paper. 

The velocity dominated matter equations can be solved exactly to give 

(2.15) °<f>(t,x) = A(x)logt + B(x) 
for given functions A and B on S and 

(2.16) °n(t,x) = A 2 (x)t- 2 

(2.17) °u a (t,x) = tlogt(A(x))- 1 e a (A(x))+tB a (x) 
for given quantities A(x) and B a (x) on S. 

2.3. Statement of the main theorems. The main theorems can now be 
stated. 

Theorem 2.1. Let S be a three-dimensional analytic manifold and let 
(° 9ab(t) , k a b(t) ,° (f)(t)) be a C w solution of the velocity dominated Einstein- 
scalar field equations on S x (0, oo) such thatttr°k = —1 and each eigenvalue 
A of —t°k a b is positive. Then there exists an open neighbourhood U of Sx {0} 
in S x [0, oo) and a unique C w solution (g a b(t), k a b(t), 4>(t)) of the Einstein- 
scalar field equations on U n (5 x (0, oo)) such that for each compact subset 
K C S there are positive real numbers (, (3, a a b , with Q < j3 < a a b , for which 
the following estimates hold uniformly on K: 

1. °9 ac 9cb = 6% + o(t a \) 

2. k\ = °k a b + o(t- l+aa b) 

3. = V + o(^) 

4. d t 4> = d t °(f) + o(t- 1+ P) 

5. V c e/(<7cb) = o(t«VC) 

6. e W) = e a (V) + o(^- C ) 

Note that the condition on tr°fc can always be arranged by means of a 
time translation and that the condition on the eigenvalues of °k a b is satisfied 
provided it holds for a single value of t > 0. The positivity condition on 
the eigenvalues together with the velocity dominated Hamiltonian constraint 
imply that A 2 must be strictly positive in the velocity dominated solution. 
Thus vacuum solutions are ruled out by the hypotheses of this theorem. If 
an analogous analysis were done for the Einstein equations coupled to other 
matter models, for instance a perfect fluid with equation of state p = kp, 
k < 1, then in many cases, including that of the fluid just mentioned, the 
matter would make no contribution to the velocity dominated Hamiltonian 
constraint, so that it would not be possible to prove an analogous theorem. 
This reflects the fact that for those matter models an oscillatory approach 
to the singularity is predicted by the BKL analysis. The scalar field is an 
exception, as is the stiff fluid which will be discussed next. 
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Theorem 2.2. Let S be a three-dimensional analytic manifold and let 
(°9ab(t)> k a b(t), iJ>(t),°u a ) be a C w solution of the velocity dominated Einstein- 
stiff fluid equations on Sx (0, oo) such that ttr°k = — 1 and each eigenvalue X 
of —t°k a b is positive. Then there exists an open neighbourhood U of S x {0} 
in S x [0, oo) and and a unique C w solution (g a b(t) , k a b(t) , /i(t) , u a (t)) of 
the Einstein- stiff fluid equations on U D (S x (0, oo)) such that for each 
compact subset K C S there are positive real numbers (, (3\, /?2, a a b , with 
C < 02 < Pi < a a b , for which the following estimates hold uniformly on K: 

1. °9 ac 9cb = S a b + o(t aa b ) 

2. k a b = °k a b + o(t- 1+aa b) 

3. n = V + o(t" 2+/31 ) 

4. u a = °u a + o{t 1+ fc) 

5. °g ac e f (g cb ) =o(t* a b -<) 

The interest of these theorems depends very much on what information is 
available on constructing solutions of the velocity dominated system. Sup- 
pose that a solution of the velocity dominated constraints is given for some 
t = to > 0. The velocity dominated evolution equations constitute a system 
of ordinary differential equations which can be solved with these initial data. 
It follows from the remarks above that the solution exists globally on the 
interval (0,oo). If we define 

(2.18) °C = -°k ab °k ab + (tr°A;) 2 - 16vr p 

(2.19) °C b = V a (°k ab ) - e b (ti°k) - 8vr°j fe 
then the velocity dominated evolution equations imply that: 

(2.20) d t °C + 2t- 1 (°C) = 

(2.21) d t °C a + t- 1 ( C a ) = ~e a (°C) 

To prove this it is necessary to use the following equations for the mat- 
ter quantities, which can be derived from the velocity dominated matter 
equations in both the scalar field and stiff fluid cases. 

(2.22) d t °p = 2(tr°fc)V 

(2.23) d t °j a = (tr° k)°j a -e a (°p) 

Since °C vanishes at t = to the evolution equation for °C implies that it 
vanishes everywhere. Then the evolution equation for °C a , together with 
the fact that it vanishes for t = t$ implies that °C a vanishes everywhere. To 
sum up, if the velocity dominated constraints are satisfied at some time and 
the velocity dominated evolution equations are satisfied everywhere then 
the velocity dominated constraints are satisfied everywhere. Thus in order 
to have a parametrization of the general solution of the velocity dominated 
system, it is enough to obtain a parametrization of solutions of the velocity 
dominated constraints. The latter question will be treated in section [7|. 
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3. Framework of the proofs 



In this section the proofs of Theorems |2.1| and |2.2| are outlined. Only the 
general logical stucture of the proof is explained here and the hard technical 
parts of the argument are left to later sections. These results will be referred 
to as required in this section. 

The first step is to make a suitable ansatz for the desired solution. This 
essentially means giving names to the remainder terms occurring in the 
statements of the main theorems. Assume that a velocity dominated solution 
is given as in those statements. Then a solution is sought in the form: 



(3.1a) g ab = °g ab + <7 O c* a V& 

(3.1b) k ab = g ac (°k c b + t~ 1+aC b K c b ) 

In the following the summation convention applies only to repeated tensor 
indices and not to non-tensorial quantities like a a b . Thus in the above equa- 
tions there is a summation on the index c but none on the index b. Matter 
fields are sought in the form: 

(3.2) <(> = <l> + tPil> 

and 

(3.3a) /i = V + ^ 2+ft ^ 

(3.3b) u a = V + t 1+ ^v a 



respectively. The Einstein-scalar field equations (2^), fl2.3| ), fl2,5| ) and ( |2.6| ) 



can be rewritten as equations for r y a bl K a b and ip. Similarly the Einstein-stiff 
fluid equations can be written as equations for j a b , K a b , v and v a . This sys- 
tem of equations (for either choice of matter model) will be called the first 
reduced system. Since 7 a fe and K a b are mixed tensors, there is no direct 
way to express the fact that they originated from symmetric tensors. Instead 
it must be shown that when the first reduced system is solved with suitable 
asymptotic conditions as t — > then the quantities g ab and k ab defined by 
the above equations are in fact symmetric as a consequence of the differ- 
ential equations and the initial conditions. When allowing non-symmetric 
tensors g ab and k ab , we need to establish some conventions in order to make 
the definition of the first reduced system unambiguous. Firstly, define g bc as 
the unique tensor which satifies g a b9 bc = 8 a c . Next, use the convention that 
indices on tensors are lowered by contraction with the second index of g ab 
and raised with the first index of g bc . This maintains the usual properties of 
index manipulations in the case of symmetric g ab as far as possible. The co- 
variant derivatives in the equations are expressed in terms of the connection 
coefficients in the frame {e a } in an unambiguous way. The definition of the 
connection coefficients is extended to the case of a non-symmetric tensor g ab 
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by fixing the order of indices according to 

(3.4) g c d^ib = \{e a {gbc) + e b (g ac ) - e c (g ab ) + i d ah g c d - lic9bd - itcSad) 

where = 6 c ([e a , e&]) are the structure functions of the frame. 

Finally, in order to define the Ricci tensor in the evolution equation for 
k a b we define R a b to be the Ricci tensor of the symmetric part s g a b = 



(l/2)(g a b +gba) of g a b- I n Lemma 5.2 it is shown that given a solution of the 



velocity dominated system as in the statement of one of the main theorems 
|2.1| -2.2, any solution of the first reduced system which satisfies points lL|] 



of Theorem 2.1, in the scalar field case, and points of Theorem 2.2, in 



the stiff fluid case, gives rise to symmetric tensors g a b and k a b, and thus to 
a solution of the Einstein- matter evolution equations. It then follows from 
Lemma [O] and the remarks following it that the Einstein-matter constraints 
are also satisfied. It follows that to prove the main theorems it is enough to 
prove the existence and uniqueness of solutions of the first reduced system 
of the form given in the main theorems. 

The existence theorem for the first reduced system will be proved using 
the theory of Fuchsian systems. Since the form of this theory we use in the 
following concerns a system of first order equations it is not immediately 
applicable, due to the occurrence of second order derivatives of the metric 
and scalar field. It is necessary to introduce some suitable new variables 
representing spatial derivatives of the basic variables. Define X bc = t^e c (j a b ) 
and, in the scalar field case, to a = t^e a (tp) and \ = tdtip + Pifii where £ 
and (5 are positive constants. The evolution equations satisfied by these 
quantities are given explicitly in ( 5. 13b| ) and ( 5.15| ). With the help of the 



new variables these equations together with the first reduced system can 
be written as a first order system. Call the result the second reduced 
system. It is easy to show, using the evolution equations for differences like 
7^ c — ^e c (7 a b ) which follow from the second reduced system, that the first 
and second reduced systems give rise to the same sets of solutions under the 
assumptions of the main theorems, together with corresponding assumptions 
on the new variables. Thus it suffices to solve the second reduced system, 
which is of first order. 

If it can be shown that the second reduced system is Fuchsian, then the 



main theorems follow from Theorem [L2|. In fact it is enough to show that 



the restriction of the system to a neighbourhood of an arbitrary point of S is 
Fuchsian. For the local solutions thus obtained can be pieced together to get 
a global solution. Moreover asymptotic estimates as in the statements of the 
main theorems follow from corresponding local statements, since a compact 
subset of S can be covered by finitely many of the local neighbourhoods. 

In sections || and ^ it is shown that the second reduced system is Fuch- 
sian on some neighbourhood of each point of S for a suitable choice of the 
constants a a b , (3, f3\, fti and C depending on the given velocity dominated 
solution. This requires a detailed analysis of the degree of singularity of 
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all terms in the second reduced system, in particular that of the Ricci ten- 
sor. The result of the latter is Lemma [T^. This is the hardest part of the 
proof. Due the above considerations, it is clear that the main theorems 



follow directly from Theorem 4.2 



4. FUCHSIAN SYSTEMS 



The proofs of Theorems 2.1 and 2.2 rely on a result of Kichenassamy 



and Rendall|13] on Fuchsian systems which uses a method going back to 



Baouendi and Goulaouic|2|]. The result of |13[] will now be recalled. It 
concerns a system of the form: 

du 

(4.1) t— + A(x)u = f(t,x,u,u x ) 

Here u(t, x) is a function on an open subset of R x R n with values in R fc and 
A(x) is a matrix- valued function. The derivatives of u with respect to the 
x variables are denoted by u x . The function / is defined on (0, T ] x U\ x f7 2 > 
where U\ is an open subset of R n and U2 is an open subset of H k+nk , and 
takes values in R fe . We assume that A(x) is defined on U\. 

In this and later sections it will be useful to have some terminology for 
comparing the sizes of certain expressions. 

Definition 4.1. Let F(t,x,p), G(t,x,p) be functions on (0, To] x U\ x U2, 
where U±, are open subsets of R n and R^ respectively. Then we will say 
that 

F =4 G 

if for every compact K C U\ x U2, there is a constant C such that 
\F(t,x,p)\ < C\G(t,x,p)\ for t £ (0,t ], (x,p) € K. 

In the particular case that G is just a function of t we will often use the 
familiar notation F = 0(G(t)) to replace F =4 G. The notation F = o(G(t)) 
will also be used to indicate that F/G tends to zero uniformly on compact 
subsets of U\ x U2 as t — > 0. 

In the theorem on Fuchsian systems the function / is supposed to be reg- 
ular in a sense which will now be explained. To do this we need the notion 
of a function which is continuous in t and analytic in other (complex) vari- 
ables. This means by definition that it should be a continuous function of all 
variables, that the first order partial derivatives with respect to all variables 
other than t should exist and be continuous, and that the Cauchy-Riemann 
equations should be satisfied in these variables. For further remarks on this 
concept sec |l5|]. Assume that there is an open subset U of C n+k+nk whose 
intersection with the real section is equal to XJ\ x U2 and a function / on 
(0, To] x U continuous in t and analytic in the remaining arguments whose 
restriction to (0, To] x Y{J l+k+nk is equal to /. The function / is called reg- 
ular if it has an analytic continuation / of the kind just described, and if 
there is some 6 > such that / and its first derivatives with respect to the 
arguments u and u x are 0(t e ) as t — > 0, in the sense introduced above. 
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For a matrix A with entries A a b let ||^4|| = sup{||vlx|| : |[x|| = 1} (operator 
norm) and ||^4||oo = m & x ab \A a b \ (maximum norm). Since these two norms 
are equivalent the operator norm could be replaced by the maximum norm 



in the statement of the theorem which follows and in Lemma 4.3 below. 
However in the proof of that lemma below the use of the operator norm is 
important. 

Theorem 4.2. Suppose that the function f is regular, A(x) has an analytic 
continuation to an open set U\ whose intersection with the real section is U\, 
and there is a constant C such that \\o~ A ^\\ < C for x € U\ and < a < 1. 



Then the equation (4-1) has a unique solution u defined near t = which is 
continuous in t and analytic in x and tends to zero as t — > 0. If f is analytic 
for t > then this solution is also analytic in t for t > 0. 

Remark 4.1. Under the hypotheses of the theorem spatial derivatives of any 
order of u are also o(l) as t — > 0. This follows directly from the proof of the 
theorem in [13|. 



Remark 4.2. If the coefficients of the equation depend analytically on a 
parameter and are suitably regular, then the solution depends analytically 
on the parameter. It suffices to treat the parameter as an additional spatial 
variable. 

The statement of the theorem is not identical to that given in |]l^] but 
the proof is just the same. We can write f(t,x,u,u x ) = t e g(t,x,u,u x ) for 
some bounded function g. By replacing t as time variable by t e , it can be 
assumed without loss of generality that 6 = 1. Then the iteration used in 
the proof given in [13] converges to the desired solution under the regularity 



hypothesis we have made on /. 

For the applications in this paper an extension of this result which applies 
to equations slightly more general than ( |4.1| ) will be required. These have 
the form 

du du 
(4.2) t— + A{x)u = f(t, x, u, u x ) + g(t, x, u)t— 

If we have an equation of this form where / and g are regular then an 
analogous existence and uniqueness result holds. For we can rewrite ( |4.2| ) 
in the form 

(4.3) 
du 

t— + A(x)u =[/-(/- g(t, x, u^jAix) + (J - g(t, x, u))" 1 /^, x, u, u x ) 

If we call the right hand side of this equation h(t,x,u,u x ) then it satisfies 
the conditions required of f(t,x,u,u x ) in Theorem [4.2| . In other words 
h is regular. For if g(t,x,u) is 0(t e ) then (/ — g{t,x,u))~ l is 0(1) and 
[I-{I-g{t,x,u))- x ] is 0(t e ). 
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In the following it will be necessary to verify the regularity hypothesis for 
some particular systems, and some general remarks which can be used to 
simplify this task will now be made. In these systems we always have 

rn 

(4.4) f(t, x, u, u x ) = tiFi(x, v(t, x),u, u x ) 

i=l 

where the F{ are analytic functions on an open subset V of Yl n+l+k+nk which 
includes U\ x {0} x U2, ii, . . . ,t m are some functions of t and v(t,x) is a 
given function with values in R'. The tj are continuous functions on (0,7o] 
which tend to zero as t — > as least as fast as some positive power of t and 
are analytic for t > 0. The function v(t, x) is the restriction of an analytic 
function on (0, To] x U\ to real values of its arguments. Each component of 
v tends to zero as t — » 0, uniformly on U±. These properties ensure that / 
is regular. For the functions F{ have analytic continuations to some open 
neighbourhood V of V in Q n+l+k+nk . 

In the examples we will meet the functions ti are either positive powers of 
t or positive powers of t times positive powers of log t. The role of v(t, x) is 
played by the functions tj, the components of the velocity dominated metric 
°g a b, their spatial derivatives of first and second order, and the components 
of the inverse metric multiplied by suitable powers of t so that the product 
vanishes in the limit t — ► 0. We saw in the last section that °g a b can be 
written in the form t 2K which, for each fixed t, is an entire function of K. 
Thus t 2K w is analytic on any region where K(x) is analytic. The velocity 
dominated metric and its derivatives all tend to zero uniformly on compact 
sets as t — > while the same is true of the inverse metric multiplied by a 
suitable power of t. 

Next a criterion will be given which allows the hypothesis on the matrix 
A in Theorem |4.2| to be checked in many cases. 

Lemma 4.3. Let A(x) be a kxk matrix-valued continuous function defined 
on a compact subset of R™. If there is a constant a such that, for each 
eigenvalue A of A(x) at any point of the given compact set, ReA > a then 
there is a constant C such that the estimate Ht^^l — Ct a holds for t small 
and positive and x in the compact set. 

Proof The general case can be reduced to the case a = by the following 
computation: 

(4.5) \\t A ' aI \\ < \\t A \\t~ a 

In the rest of the proof only the case a = will be considered. Without 
the parameter dependence the result could easily be proved by reducing the 
matrix to Jordan canonical form. The difficulty with a parameter is that 
the reduction to canonical form is in general not a continuous process. At 
least it can be concluded from the continuity properties of the eigenval- 
ues that there is a (3 > such that all eigenvalues satisfy ReA > f3. Let 
s = log t. Then the problem is to show that for a fixed sq there is a constant 
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C > such that ||e _Sj4 || < C for all s > sq. By scaling t we may suppose 
without loss of generality that sq = 0. For each x we can conclude by reduc- 
tion to canonical form that exists a value s x of s such that the inequality 
|| e — a a A(x)|| ^ e -/3s x /2 fooi^g gy continuity of the exponential function there 
is an open neighbourhood U x of x where this continues to hold for the given 
value of s x . Let C x = sup{||e _,sA ^ || : s G [0, s x ],y G U x }. It follows that for 
any s G [0, oo) and any y G U x we have 

(4.6) ||e~ sA(?/) || < C x \\e~ SxA(y) \\ [s/Sx] < C x e' 0Sx[s/sx]/2 < C x 

By compactness, it is possible to pass to a subcover consisting of a finite 
number of the sets U x and letting C be the maximum of the corresponding 
C x we obtain the required estimate. 

Remark 4.3. More generally, an analogous estimate is obtained if A(x) is 
the direct sum of a matrix B{x) whose eigenvalues have positive real parts 
and the zero matrix. This is obvious since in that case t A ^ is the direct 
sum of t B<yX ^ and the identity. 

5. Setting up the reduced equations 

In this section we introduce the adapted frame {e a } and the auxiliary 
exponents {q a } which will be used in the curvature estimate. 

Let xq be given. Let °g a b, °k a h , be solutions of the velocity dominated 



6' 

evolution and constraint equations ( [2.1l| ) and ( 2.1C| ). Let p a be the eigenval- 
ues of K\ = -t°k a b . Assume that {p a } are such that p a (xo) > 0, a = 1, 2, 3, 
^2 a p a = 1 (Kasner condition) and p a are ordered so that p a < pt, for a < b. 
Fix an initial time to G (0, 1). We will in the following restrict our consider- 
ations to t G (0, 1). 

If K has a double eigenvalue at xq, then in general the eigenvalues and 
eigenvectors of K are not analytic in a neighbourhood of xq, and therefore in 
general it is not possible to introduce an analytic frame diagonalizing K in 
a given neighbourhood. We will avoid this problem by using the well known 
fact that if p a > is a double eigenvalue of K(xq), the eigenspace of the pair of 
eigenvalues corresponding to p a > is analytic in a neighbourhood of xq. This 
means in particular that it is possible to choose an analytic frame which is 
adapted to the eigenspace of the pair of eigenvalues. This will play a central 
role in what follows. 

Choose numbers ao, e > so that e = ao/4 < min{p a (xo)}/40. 
1. Cases I, II, III: We will distinguish between the following cases: 
(I) (near Friedmann) maxa^ \p a — pi,\ < e/2, a = 1, 2, 3. 
(II) (near double eigenvalue) max a 6 \p a — Pb\ > e/2, and \p> a — Pb>\ < 
e/2 for some pair a', b' , a' ^ b'. Denote by p± the distinguished 
exponent not equal to Pa'iW- 
(III) (diagonalizable) min^t, \p a — pb\ > e/2 

By reducing e if necessary we can make sure that condition I, II or III 
holds at xq if the maximum multiplicity of an eigenvalue at xq is three, 
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two or one respectively. The conditions I, II, III are open, and hence 
there is an open neighbourhood Uq 3 xq such that, if condition I, II 
or III holds at xq, the condition holds in Uq and further, for x € Uq, 
mm. a {p a (x)} > 20e. 

2. Auxiliary exponents {q a }: Let Uq 3 xq be as in point |l|. We will 
define analytic functions q a ,a = 1,2,3, called auxiliary exponents, in 
Uq, with the properties 

(a) q a > (positivity) 

(b) q a < % if o, < b (ordering) 

( c ) Ea Qa = 1 (Kasner) 

The auxiliary exponents q a will be defined in terms of the eigenvalues 
p a of K depending on whether at xq we are in case I, II, or III. 

Let q a {xo) = p a {xo), a = 1,2,3, and choose q a on Uq satisfying the 
positivity, ordering and Kasner conditions such that in cases I, II, III, 
the following holds. 
(I) 9a = 1/3, a = 1,2, 3 
(II) q± = p±, q a > = q b > = \{l - qj_) 
(III) q a = p a 

Then q a are analytic on Uq. Note that it follows from the definition of 
q a that qi > minjjpi} and that max a \q a — p a \ < e/2. 

3. The frame {e a }' In each case I, II, III define an analytic frame {e a } 
with dual frame {9 a }, by the following prescription: {e a } is an ON 
frame w.r.t. g(to), and in case II, III the following additional condi- 
tions hold. 

(II) ej_ is the eigenvector of K corresponding to q± and e a > , ey span 

the eigenspace of K corresponding to the eigenvalues p a ',Pb'- 
(III) e a are eigenvectors of K corresponding to the eigenvalues q a - 

We will call {q a }, {^a}, {^ a }j satisfying the above conditions adapted. In 
the following we will work in a neighbourhood Uq defined as above and 
assume that we are given adapted {q a }, {e a }> {# a }> on Uq. 

The role of a a b will be to shift the spectrum of the system matrix to be 
positive. We need to choose ao larger than e to compensate for the fact that 
the q a are not the exact eigenvalues of K. 

In the following T a b will denote frame components of the tensor T w.r.t. 
the frame {e a }- Define the rescaled frame e a = t~ qa e a with dual frame 
Qa _ fq a Qa an( j ,-[ eno t e by j 1 ^ trig e a frame components of the tensor T. 
Then we have 

(5.1) fab = t-^Tab, f% = t"n. 

It follows from the definitions, that in case III, 



(5.2) 



K% = 5\q b , °k a b = -r X 5 a b q h 
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while in case II, the tensors K, °k and °g are block diagonal in the frame 

M, 

(5.3) K a ' x = 0, °k a ' x = 0, V± = 
For s € R, let (s) + = max(s,0), for a, b G {1,2,3}, let 

(5.4) a a b = 2(q b - q a ) + + a , 

and a a b = \q b — q a \ + (Xq. In view of the relation 2(s) + — s = \s\ we have 

(5.5) a\ + q a -q h = a\. 



The following identities are an immediate consequence of equations ( |5.2| ) 
and (|5.3|), together with the fact that in case I, = 1/3, a = 1, 2, 3. 



(5.6a) °g ab t ab c = °g ab t^ 

(5.6b) °g ab ti» = °g ab t^ 



Note that in ([5l]) no summation over indices is implied. 

Lemma [4.3| implies the following estimate for the rescaled frame compo- 
nents of °g. 

Lemma 5.1. 

(5.7) \\°~9ab\\oo<Ct- e , HVIIoc <Ci~ e 

Proof. Let Q b = 5 a b q b . A direct computation starting from the matrix form 
of the velocity dominated evolution equation gives 

td t G ab = 2G ac {K c b - Q c b ) 

which has the solution 



G a b(t) — G ac (to) 



to 



From the definition of q a and the frame e a the spectrum of K — Q is of 
the form p a — q a . Therefore since \p a — q a \ < e/2, we find that the spectrum 
of 2(K — Q) is contained in the interval (— e,e). Therefore it follows from 
Lemma (D| that \\G\loo < Ct~ e . 

Similarly, using the fact that G _1 satisfies the equation 

td t G-J = 2(Q a c - K a c )G~ b 

an application of Lemma 4.2 yields the estimate | IG^ - 1 1 |oo < Ct~ e . □ 

We are now ready to describe the ansatz which will be used to write 
the Einstein-scalar field and Einstein-stiff fluid systems in Fuchsian form. 
Assume that a solution of the velocity dominated constraint and evolution 
equations, °g a b, °k ab , with adapted frame, coframe, and auxiliary exponents 
{e a }, {6 a }, {q a }, is given. Let a a b be defined by 
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Let C = e/200. We will consider metrics and second fundamental forms 
g, k of the form 

(5.8a) g ab = °g ab + V c i Q V b 1% = o(l) 

(5.8b) g ab = °g ab + t aa ^ a c °g cb j a c = o(l) 

(5.8c) e c ( 7 a fe ) = i"^ A£ c = o(l) 

(5.8d) fc a6 = <? ac (°£; c b + t- 1+QC ^ c fe ) K c fe = o(l) 

The form (5.8b) is a consequence of ( |5,8a| ). To see this, note that °g ac g c b = 
5 a b + t aa b~j a b and that g ac0 g c b = ^g ac g c b)~ l - Thus the desired result follows 
from the matrix identity (7 + A)" 1 = I — A + (I + A)^ 1 A 2 and the fact that, 
using 2{x) + — x = \x\, it can be concluded that 

(5.9) a a e + a e b - a a b = \q e - q a \ + \q b - q e \ - \q b - q a \ + a Q > a . 

The latter relation shows that each component of the square of 7% vanishes 
faster than the corresponding component of 7°;, itself. 

Let f3 = e/100. In addition to ( |5.8j ) we will use the following ansatz for 
the scalar field 

(5.10a) c/> = V + t /3 V ip = o(l) 

(5.10b) e a (i/>) = t'^oJa uj a = o(l) 

(5.10c) tdfif; + Pip = x X = o(l) 

and for the stiff fluid case, 

(5.11a) fi = °fi + t- 2+l3l u i/ = o(l) 

(5.11b) u a = °u a + t 1+ ^v a v a = o(l) 

Equations |[|), ( |57TO| ) and (|5TlD with the exception of ( |5T8cD and ( |5T0bD 
will be used to derive the first reduced form of the field equations, and 
equations ( |5.8c| ) and ( 5.10b ) for the spatial derivatives of 7 a fe and ip, will be 
used to derive the second reduced system. 

Note that in view of (5.1), we have 

(5.12a) 9ab = 9ab + °9act &C ^ c b 

(5.12b) g ab = °g ab + t &a ^ a c °g cb 

(5.12c) ~k ab = ~gac(°~k c b + t- 1+&c *K c b ) 

We use the following conventions throughout: 

• indices on velocity dominated fields °g a bi k ab ^u a are raised and low- 
ered with °g ab , while indices on other tensors are raised and lowered 
with g ab . 

• the dynamic tensor fields J a b ,K a b in g ab ,k ab are always used in mixed 
form and only in {e a } frame components. 

• the dynamic 1-form v a in the velocity field u a is always used with lower 
index. 
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5.1. The reduced Einstein matter system. In this section, we describe 
the first reduced system for the Einstein-scalar field evolution equations, 
derived from ( |2.3j ) using the ansatz given by equations ( |5.8[ ) and (5.10) 
for g a b,k a b,(p in terms of the velocity dominated solution g a bi kab, 4>, the 
auxiliary exponents {q a } and the dynamical fields j a b , K a b , ip. Similarly we 
describe the first reduced system for the Einstein-stiff fluid evolution equa- 
tions obtained using the equations (5.11). For convenience we use the term 
'Einstein-matter system' to describe the Einstein-scalar field and Einstein- 
stiff fluid system collectively. 

The tensor g ab of the form ( |5,8a| ) is not a priori symmetric, but it will fol- 
low from Lemma 5.2 that the solution to the Fuchsian form of the Einstein- 



matter evolution equations will be symmetric. It is convenient to introduce 
the symmetrized tensor 

1 



g a b 



:(g a b + gba)- 



Let R a b be the Ricci tensor computed w.r.t. the symmetrized metric g a b, 
see section || for details. 

By substituting into the evolution equations (|2.3|) with R ab replaced by 
R a b, defined in terms of "y a b and X a bc , we get the following system for 



s 



ryO- K a \a 
7 fe) K b) A be- 

(5.13a) 



td a a b + a a bl a b + 2k\ + 2 7 a e (t fc e 6 
(5.13b) 

tdAZ^t^tda^ + tfe^b) 

(5.13c) 

td t K\ + a a b K a b - (t°k a b )(tvn) = t ao {tiK)K a b 
(5.13d) +t 2 - aa b( s R a b - M%) 
where M ab is is given by 
(5.14a) M ab = 87re a (0)e 6 (0) 
(5.14b) M ab = IQ-Kjiu a u h 



nt°k a e h e b 



-2t aa e +a£ »- 



b l e« b 



for the Einstein-scalar field system 
for the Einstein-stiff fluid system 



M a b will be estimated in sectio n |6| N ote that the power of t occurring on 
the right hand side of equation (5.13a) is positive due to (5.9). 

The wave equation ( |2 .6| ) becomes the following system of equations for 

(5.15a) td t if> + f3ip - x = 

(5.15b) td t u a = t C [e a ( X ) + (C - PKW] 

(5.15c) td tX + (3x = t a °-PtiK(A + t p x) + i 2_/3 A°0 + t 2 ^V a uj a 

Let U = (7 a b , K a b , ip, X-, A^ c , uj a ). Then we can write the second reduced 
system in the Einstein-scalar field case, which consists of equations ( |5.13 ) 
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and ( 5.1 5| ) in the form 

td t U + AU = T{t,xMMx)- 

for a matrix A and a function T. We will prove that this system is in 
Fuchsian form. 

The Einstein-stiff fluid equations will be treated in a similar way. How- 
ever, due to the complexity of the equations in that case, they will not be 
written out more explicitly than is absolutely necessary to understand the 
essential features of their structure. The second reduced system in the stiff 
fluid case can be brought into the generalized Fuchsian form 

(5.16) td t U + AU = F(t, xMMx) + Q(t, x,U)d t U 

already introduced in section |j. In order to do this it is useful to introduce 
some abbreviations for certain terms in ( |2.9[ ) so that the equations become 

(5.17a) %t - 2(trk)fi = -2\u\ 2 d t fi - \[iu a d t u a + F x 

(5.17b) 8 t u a + (trk)u a + (l/2) / u- 1 e a (^) = - 2(trk)fj]u a 

+ (1 + M 2 )- Vu a d iUfe - [(1 + H 2 )- 1 ' 2 - l]/i _1 ea(M) + F 2 

The expressions F\ and F 2 contain only terms which can be incorporated 
into T in (5.16). Next the ansatz ( |5,11| ) must be substituted into these 



equations. The result is: 
(5.18a) 

td t u + fiiv = -2t 3_/?1 (1 + ttxkfn + 2(1 + tlrk)v + t 3 '? 1 [d t fi - 2(trk)/j] 
(5.18b) 

td t v a + (3 2 v a = -(l + ttik)v a + t~^[d t u a + {trk)u a + (l/2)V _1 e a (V)] 

The expressions on the left hand side of the above form of the Euler equations 
written in terms of the basic variables \jl and u a occur on the right hand 
sides of the above evolution equations for v and v a . In order to get a fully 
explicit form it would be necessary to substitute for these expressions and 
then express the final result in terms of v and v a . This is, however, neither 
necessary nor even helpful for the analysis to be done here. 

Next we will consider the matrix A and prove that A is a direct sum of 
a matrix with spectrum bounded from below by a positive number, with a 
zero matrix. (The arguments in the scalar field and stiff fluid cases are very 
similar.) It is therefore of a form such that the theory presented in section 
§ applies. In addition we must show that !F(t,x,U,U x ) = 0(t s ) for some 
5 > and, in the stiff fluid case, that Q(t,x,U) satisfies a similar estimate. 
This will be done in the next section. 

The matrix A is block diagonal and therefore it is enough to consider each 
block separately. The rows and columns of A corresponding to X a bc , uj a are 
zero, and therefore this A is the direct sum of a matrix corresponding to 
7, k, if), X; with a zero matrix in the scalar field case and the direct sum of a 
matrix corresponding to 7, k, u, v a with a zero matrix in the stiff fluid case. 
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We now consider the spectrum of this matrix. The submatrix corresponding 
to 7, k is upper block triangular. The 7, 7 block is given by 

1 a b ^a a bl a b + 2[ 1 ,t°k] a b 

To estimate the spectrum of this, it is necessary to consider the cases I, II, III 
separately. Working in a frame which diagonalizes k, t°k a b = —5 a b pb, and 
hence in this case 

2[ 1 ,t°k] a b = -2(p b -p a ) 1 a b 
Therefore, in case III, we get using the definition of a a b , 

a a bl % + 2[7, t°k] a b = 2((p b - p a ) + - (p b - p a ) + a )7 a fo 

and hence using — x > for all x € R, the spectrum of the 7,7 block 
is bounded from below by ao in case III. 

Next consider case I. In this case, a a b = ao = 4e and the spectrum of 
7 a fe 1— > 2(7, t°k] a b is bounded from below by e, which shows that the spectrum 
of the 7, 7 block is bounded from below by 3e in case II. 

Finally, in case II, °k a b is block diagonal in the adapted frame. The 
spectrum of "y a b 1— > 2[y, t°k] a b consists of 2(p a i — py), 2(p a ' — P±), and 0. 
Now using the definition of a a b for case II and arguing as above, we get the 
lower bound 3e for the spectrum of the 7, 7 block in case II. Therefore the 
spectrum of the 7, 7 of A is bounded from below by 3e. 

Next we consider the k, k block. This is of the form 

K% I ► a a b K a b - {t°k a b )tiK 

First consider the action on the trace-free part of K a b . Then the spectrum 
is given by a a b > ao. On the other hand, restricting to the trace part of 
n a b , which is diagonal, we see that the spectrum is ao + 1. Therefore the 
spectrum of the n, k block is bounded from below by ao- 
The ip, x block is of the form 



which has spectrum f3 > 0. The v, v a block is diagonal with eigenvalues f3± 
and fa- 

Therefore, in view of the facts that 3e > (3 > 0, f3\ > and $2 > 0, the 
desired properties of the spectrum of A have been verified. 

Given a solution U = (7 a fe , n a b , tp,x> ^tc w a) °f * ne reduced system for the 
Einstein-scalar field equations, define g ab , k ab and 4> by (|5.8|) and ( |5.10D , 
Similarly, given a solution U = (^ a b , K a b , v a X bc ) of the reduced system for 
the Einstein-stiff fluid equations, define {g a b,k ab , fi,u a ) by ( |5.8| ) and ( |5.11 ), 

If it can be shown that g ab and k ab are symmetric then a solution of the 
Einstein-scalar field equations is obtained. The next lemma gives sufficient 
conditions for this to be true. 

Lemma 5.2. Let a solution of the velocity dominated Einstein-matter sys- 
tem be given on S x (0, 00) with all eigenvalues of —t°k a b positive and 
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ttr°k = — 1. Let U be a solution of the reduced system for the Einstein- 
scalar field or Einstein-stiff fluid system corresponding to the given velocity 



dominated solution, with U = o(l). Define (g a bi k a b, (/)) by (5.6) and (5.1C) 



in the scalar field case, and define (g a b, k ab , L 1 -, u a) by (5.1 ) and (5.11) in the 
stiff fluid case. Then g ab and k ab are symmetric. 

Proof. From the evolution equation for 7 a b and the definitions of g ab and 
k ab it follows that dtg a b = —2k ab . Similarly an equation close to the usual 
evolution equation for k a b can be recovered from (|5.13d ). It differs from the 



usual one only in the fact that R a b is replaced by s R a b - From these equations 
we can derive the equations: 

(5.19a) d t (g a b - 9ba) = -2(k ab - ha) 

(5.19b) d t (k ab - k ba ) = (tnk)(k ab - k ba ) 

It follows from the assumptions on ^ a b and n a b together with the definition 
of k ab , that the components of k ab — k ba are o(t~ 1+v ) for some r/ > 0. Hence 
the quantity il ab = t ^(kgjt, — k ba ) tends to zero as t — > 0. It satisfies the 
equation: 

(5.20) td t n ab + 7]Q ab = (ttvk + 1)0^ 



From Theorem 4.2 we conclude that Q ab = 0. Thus k ab is symmetric. It 



then follows immediately from (5.19a) and the fact that g ab = o(l) that g ab 



is also symmetric. □ 

6. Curvature estimates 

Let s R a b be the Ricci tensor computed w.r.t. the symmetrized metric 
S 9ab = \{§ab + 9ba) ■ I n order to get a Fuchsian form for the Einstein- 
matter evolution equations, we need the following estimate for the frame 
components of S R a b , 

(6.1) t 2 ~ aa b s R a b = 0(t 5 ), for some 5 G (0, e). 

In doing the estimates we will use the notion of comparing the size of 



functions introduced in Definition 4T . In proving that the second reduced 
system 

tdtu + A(x)u = /(i, x, u, u x ) 
is in Fuchsian form, one essential step is to prove an estimate of the form 

f4t S 

for some 5 > 0. In the present section, we accomplish this task for the expres- 
sion t 2 ~ aa b s R a b , which is now considered as a function r(t,x,v(x,t),u,u x ), 
where v(x,t) is defined in terms of the solution °g a bi°k ab to the velocity 
dominated system and the data {e a }, {9 a }, {q a }, etc. defined in section [|, 
and u consists of the variables 7 a b , X a bc . In terms of the relation =<I, the goal 
is to prove 

(6.2) t 2 ~ aa b s R a b 4 t 5 , for some 5 E (0, e). 
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By assumption, «o = 4e, so a a b — 2e = 2(q b — q a )+ + ao/2- Therefore, the 
arguments that apply to a a b also apply to a a b — 2e. 
The symmetrized metric tensor satisfies 



0~ , + a a .-2e0~ S„,c 
gab = 9ab + t b 9ac 7 6) 



(6.3a) 
(6.3b) 

To see this, note the identity 



Vfc = o(l), 



I c 



0(1) 



1 



which in view of Lemma 5.1 shows that J c b = o(l). The argument that 



S»,a 



7" c = o(l) is the same as for 7 a c = o(l). 
It is convenient to estimate the rescaled frame components. Note R a b 
t~ qa+qt > R a b . Hence in view of ( |5.5[) we need to consider 

±2-a a , r>a 



Using Lemma |5.l| and (|5.8| ) gives 

\\t 2 - &a »R a b \\oo<Ct-'\\t 2 - &a b R ab{loo 
To see this, we compute using ( |5.6a ) and ( 5.12b|) 

1 1+2— a a , £>a n 1 1+2— a a , ~ac £> n 

1 1* b R felloo = \\t b g Rcb\\oo 

ll+2-a a ,0 ~ac 5 n . 1 1 + 2-a a . +a a „ -a ~dc 5 i 

< \\t b g Rcb\\oo + \\t b id 9 R m 



<C\\t 



2-e-a' 



>>R. 



cb 1 1 oo 



where we used the triangle inequality in the form — a a b + a a c > — a c b . 

Let = # c ([e a , e&]) be the structure coefficients of the frame {e a }. The 
structure coefficients ^ c ab = 9 c ([e a , e b ]) of the frame e a are given by 



(6.4) 



Tab 



It is convenient to define T abc 
of g a b by 



^~"7a6 " log(t)(t-^e a ( % )«5 fc c - t-^e b {q a )5l) 

(Vg a e 6 ,e c ). Then f afec is given in terms 



(6.5) 



2f a6c = e a (# bc ) + e 6 (<7 ac ) - e c (^ a& ) 

+ Tiffed - 7ac5ta - 7bc5ad 



and Rdcab is given in terms of f abc and 7 a bc by 

Rdcab = e a tbcd ~ e&f 'acd — 7 ab^fcd 

~ g f9 ^bcf^adg + g f9 ^acf^bdg 
(6-6) = (Rl)dcab — (Rtydcab ~ (Rtydcab ~ (Rtydcab + (Rtydcab 

Let 

"0, ifa = 6 



1, if a / 6 
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Define 

Zabc(t) = t~*> + + t-*° 

Note that {Rl)dcab, ■ ■ ■ , (R^>)dcab, z ab-> Z abc are not tensors. In the rest of this 
section, the frame {e a } is fixed and the estimates will be done for tensor 
components in this frame. 

We will use the following lemma as the starting point of the estimates in 



this section. 








Lemma 6.1. 








(6.8a) 


9ab 




-e 




(6.8b) 


~ab 




-e 




(6.8c) 


e-afjbc 


=4 f 


-9a-2e 




(6.8d) 


e a h9cd 


4 f 


-Qa-Qb- 


3e 


(6.8e) 


~ g ab t ~q b 


4 f 


-e—q a 




(6.8f) 


e c {t qa 9ab) 


4 t~ 


-Qc+qb- 


2e 


(6.8g) 


g ab z abc 


4 t~ 


-e(j.-q h 


+ t~ qc + t qc ~ 2qh 


(6.8h) 


Tab~9 M 


4 r 


-2e( t q c - 


-qa-q d + t ~q a + f 



h = min(a, b) 



-qd~] 



Proof. The inequalities ( |6.8a| ) and ( 6.8b| ) are immediate from Lemma 5.1, 
( 5.12| ) and definition 4.1. Recalling that the variables u,u x occurring in the 
second reduced system contain 7 a fe , A^ c and its first order derivatives gives 
( |6.8c ) and d6.8d| ). (Here the inequality ( < e has been used.) 

The estimate ( |6.8e[ ) follows from Lemma 5.1, (5.6) and ( p. 12 ) together 
with the triangle inequality in the form —q b < —q a + \q a — qt>\. The estimate 
( |6.8fj ) follows from (|5.6| ) and ( |5.12| ) together with the observation that since 
°g ab is block diagonal for x £ Uq, e c °g ab is also block diagonal. The estimates 
( Mjh and ( pg) follow in a similar way starting from (|6.7| ) and fl6.4j ). In 
( 6.8h|) a log(i) term is dominated by t~ e . □ 



The following lemma gives estimates of T abc in terms of Z, 



abc 



Lemma 6.2. 






(6.9a) 


r abc ^4 t 


A abc 


(6.9b) 


r a &6 ^ t 


-2e-<? a 


(6.9c) 




-2e( t -q a +t -9fc) 


(6.9d) 


r aba ^4 t 


- 2t (t~ qa +t~ qb ) 


(6.9e) 


&a^bcd ^ t 


- 3t - Qa z bcd 


Proof. First observe that j bc =4 z bc . 


From (6.4) and (|6.8|) we have 



¥ bc g f a 4 t- 2 %t"- q < Zbc + r» + t- q < 
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Now noting e a (g bc ) 4 t~ qa ~ 2e , ( |6.9a| ) follows. 
To estimate T a bb we compute 

f abb = \t- Qa e a g bb 4 t- 2e - qa 



The estimates for T aa b,T a t> a follow directly from ( |6.9a|) and the definition of 

Zabc- 

Finally we consider ( |6.9e ). Expanding out e a Tb c d, we see that it contains 
(up to permutations of the indices) terms of the form e a eb(g c d), (^aj^dfb 
and rf cc fia>9fb- The first and second type of terms are estimated using ( |6.8d| ) 
and (|6.8ej), using the form of 7^. Finally, the third type of term is estimated 
using ( |6.8f| ) . □ 

An important consequence of the Kasner relation ^ a q a = 1, is 

(6.10) 2 + 2(gi -q 2 - q 3 ) = % 
which implies 

(6.11) t W{Qj-Qk-qi) ^ fox at i east one f / is different from 3, 

The strategy will be to eliminate as much as possible the occurence of re- 
peated negative exponents, in order to be able to use this relation. 
We make note of the following useful relations. 

(6.12a) Z abc 4 t^-" 

(6.12b) Z hhc 4 t- qh + ^ 

(6.12c) e a Z bcd 4 t~ e Z bc d 

The estimates (|6.12a ) and ( 6.12b ) are immediate from (|6.?1 ). 

For the rest of this section, we will assume that g ab is symmetric and is of 
the form given by ( 5.8a ). Let R ab be the Ricci tensor defined with respect 
to g ab . Under these assumptions we will estimate R a b . The estimate then 
applies after a small modification to S R a b - 

We now proceed to estimate the rescaled components of the Ricci tensor 
Rad = 9 bc Rdcab- Corresponding to the terms (-R1), . . . , (-R5) we have 

R a d = (Ricl)ad - (Ric2) arf - (Ric3) arf - (Ric4) ad + (Ric5) ad 

We make the following simplifying observations. 

• By the symmetry R a d = R^a we can assume without loss of generality 
that a < d. 

• Rdcab is skew symmetric in the first and second pair of indices, therefore 
we can assume without loss of generality that c 7^ d, b ^ a. 

Therefore, in the following, we can without loss of generality use the 
following convention: The indices a, b, c, d satisfy the relations 

(6.13) a < d, c 7^ d, a 7^ b 
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We will now estimate R a d by considering each term (Ricl) at ;, . 
in turn. The estimate we will actually prove is of the form 

which will imply the needed estimate for s R a d - 
6.1. (Ricl). 

(Ricl) arf = g bc e a t bcd 

where we are summing over repeated indices. Let F = t 2+Qa ~ qd (Ricl) a d 
Using Lemma and ( 6.8gj ), we have 

f 4 t 2 - u t qa - qd r qa (t- qc + t~ qd + t qd ~ 2qc ) 
4 t 2 -^(r qd - q ° + t~ 2qd + r 2qc ) 
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(Ric5) ad 



6.2. (Ric2). 



(Ric2) ad = g bc e b F, 



acd 



Let F = t 2+qa - qd (Ric2) ad . By Lemma |J and ( gT8p 5 

F 4 t 2 - A H q «- qd - q <Z acd 
By ( |6.13|) , d ^ c and using (|6.12a| ) and ( 6.11| ) this gives 

F 4 t 49l " 4e ^ t 4qi ~ 6e 
which is the required estimate. 

6.3. (Ric3). 

(Ric3) od = g^f fcd 
Let F = t 2+qa - qd (mc?,) ad . We estimate using ( |6.8hD and Lemma |6]2] 

F 4 t 2 - ie t qa - qd {t q f- qa - qc + t~ qa + t~ qc )Z fcd 

4 t 2 -^(t q f~ qd - qc + r qd + t qa - qd ~ q c)z fcd 

use c^d by (|6Tl3|) , ($12^) and 

^ t iqi ~ ie 4 t Aqi - 6e 

6.4. (Ric4). 

(Ric4) ad = ~g bc ~g f9 t hcf t adg 
Let F = t 2+qa ~ qd (Ric4:) ad . We have by Lemma 

In case a = d this gives, with h = min(6, c), m G {/, g}, using ( |6. 12| ) and 

(ED, 

t 2 (Ric4) aa ^ t 2 ~ 6t (r qh + t^ 9 " 1 + t qm ~ 2qh ){t- qa + 

, t 4 9l -66 
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Next we consider the case a < d. In case g = d, Lemma 3.2 together with 
( |6.8c ) and (|6.8gj ) gives 



F 4 t 2+qa - qd ~g bc ~gf d t bcf r 2e - q 



4 t 2 -^(t- q f- qh + t~ 2q f + r 2qh ) 

, t 4 qi -6e 



h = min(6, c) 



In case a = g we have arguing as above 

F 4 t 2+qa ~ qd ~ g bc g }a t bcf r 2e {r qa + r qd ) 

4 t 2 - 5e g fa (t- qh + t~ q f + t q i- 2qh )(r qd + t qa ~ 2qd ), h = min(6, c) 
^ t 2 ~ 6e (t- qh + r qm + t q ™- 2qh ){r qd + t q ™~ 2qd ), m = min(o, /) 

From h = min(6, c) and c ^ d which holds by ( |6.13j ), we find that either 
h < 3 or d < 3 must hold. Using this it follows using (ETTll) that m case 
a = g < d, F =4 t 4qi ~ 6e . 

It remains to consider the case when a < d and a, d, g are distinct. In this 
case, the estimates used above give 

F 4 t 2+qa - qd ^{r qh + r q s + t q °- 2qh )Z adg , h = mxn(b, c) 

t 2 - 6e (t~ qh + *-» + t q 9- 2qh ) 

(fid _|_ fQa-2q d _|_ ^la-qd^Hg _|_ j-2q a -2q d -q g _j_ ^.~q g _|_ ^q g -2q d ^ 

By construction, /i < 3 or d < 3 must hold, which in conjunction with the 
fact that in the present case, g ^ d gives using ( 6.11 ), 

F 4 t 4<7l_6e 

The above proves that the required estimate i 2+<?a_<?d (Ric4) a ^ 
holds. 



6.5. (Ric5). 

(Ric5) ad = g bc g f9 f acf f bdg 

The estimate for (Ric5) a ^ is the most complicated, and will be done in 
several steps. We review the steps which will be used here. In each step the 
conditions on the indices a, c, /, b, d, g which leads to the required estimate 
may be excluded from our considerations. Recall that a < d may be assumed 
and also note by ( |6.13| ) we may assume without loss of generality that a ^ b, 
c + d. 

The steps we will use are: 

1. a = d can be excluded, so a < d may be assumed. 

2. g = d can be excluded, so g ^ d may be assumed. 

3. g = a can be excluded, so g ^ a may be assumed. 

4. b = d can be excluded, so b ^ d may be assumed. 

5. a = c can be excluded, soa^c may be assumed. 



QUIESCENT COSMOLOGICAL SINGULARITIES 



27 



When all the above claims are verified, we may restrict our considerations 
to the indices satisfying the conditions 

(6.14) a < d, a j^z b, c ^ d, g ^ d, g ^ a, b ^ d, a ^ c 

These conditions imply that the indices {a,d,g}, {a,d,c} and {a,d,b} are 
distinct, so ( 6.14j) implies g = b = c, as all indices take values in {1,2,3}. 
Therefore the required estimate for (Ric5) a d will hold if we can verify that 
it holds under ( |6.14 ) in conjunction with the condition g = b = c, which is 
the final step. 
Let 

F = £ 2+9a ~ 9d (Ric5) d. 



Case a = d: In case a = d, Lemma 6.2 and ( 6.12a| ) give 

p _^ ^2-6e^2(gi-g 2 -'?3) 

Therefore we may assume a < d in the following. Further by ( |6.13 ), a ^ b 
and c ^ d. 

Case g = d: Next consider the case g = d. Then using Lemma 3J2, (|6.7| ) 
and ( 6.8e ) we have 

F = t 2+ ^g bc 9 fd T acf t bdd 

4 t 2+ "~g bc ~gf d t acf t- 2 ^ 

^ j2-%e+q a -q d ~q c U~q a _|_ £-<?c _|_ £-<?d _j_ £<?a-<? c -<?d _|_ fQc-qd-Qa _|_ fld-qa-qc^ 
j.2-6e^-q d -q c _|_ j.q a -q d -2qc 

_|_ f.q a -2q d -q c _|_ j2(q a -q d -q c ) _|_ £-2<? d _j_ £-2<?c^ 

By ( |6.13| ), d c, this gives F =<! t 4< ?i— 6e ^ n cage g = d. 

Case g = a: Next consider the case g = a. Then we have 
F = t 2+ ^~g hc ~gf a t acf t hda 



use Lemma 6.2 and (6. 



^ t 2-5e t q a -q d ~bc( t ~ qa + t ~q c + t^' 2 ^) Z bda 



use Lemma 5.1 



use (|T2a|) and a < d, d / c 
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At this stage we may assume 
(6.15) a < d, a j^z b, c / d, g ^ d, g / a 



Case b = d: Next consider the case b = d. In this case 
F = t 2+ "g dc g^t acf t ddg 

4 t 2+qa - qd g dc g f9 f acf t- 2e (t- qd + r q °) 
4 $2-6^.-^-90 + t qa - qd - qs )Z acf 

use d ^ c and d ^ g from fl6.15| ) 
^ t iqi ~ 6t 



Case a = c: Next consider the case a = c. In this case we have using 
Lemma 3.2 

F ^ t 2+qa - qd g ba g f9 t- 2e {t~ qa + t- q f)t- 2e Z bdg 

use ( |6.8e| ) 

^ t 2 - 6e (t- qd + t qa - qd ~ q o)Z Wg 
use g / d from (|6~T5| ) and ([Tl2a|) 

At this stage we may assume 

(6.16) a < d, a j^z b, c ^ d, g ^ d, g ^ a, b ^ d, a ^ c 

As discussed above, if we can prove that the required estimate holds under 
the condition g = c = b we are done. 

Case g = b = c: Next consider the case g = b = c. In this case, after 
making the substitutions g = c and b = c, 

f = t 2+q «- q *g cc g fc r acf r cdc 



use Lemma 6.2 

(6.17) ^ t 2 - 3e t qa ~ qd g fc f acf (r qd + t" 9c ) 

To estimate F we must now consider the cases f = d, f = c, f = a sepa- 
rately. 

Case g = b = c and f = d: In case g = b = c and / = d, we have from 
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use ( |6.8e| ) and Lemma [T^ 

^ t 2-6 Hqa - qd - gcZacd 

use d from ( 6.16j ) and ( |6,12a| ) 

Therefore we may exclude condition ( |6. 16|) in conjunction with f = d from 
our considerations. 

Case g = b = c and / = c: In case g = b = c and / = c we have from 

o 

F 4 t 2 - 3e t qa - qd g cc f acc (t- qd + t~ qc ) 



use Lemma 6.2 



Case g = b = c and f = a: The only remaining case is 5 = 6 = c and 
f = a. In this case we get from ( |6.17D using Lemma |6.2| 

F ^ t 2 - 5e t qa - qd g ac Z aca (t- qd + t~ qc ) 



use ( |6.12b|) and 



_^ ^2— 6e/^g a — 2g<j _|_ £?a-?d- "9cV— 9a 



Therefore it now follows that under ( 6.16 ), the required estimate 

F t 49l " 6e 

holds and hence by the above argument it follows that this estimate holds 
under ( |6.13|) . 



This proves for a symmetric metric satisfying (5.8) the estimate 
(6.18) t 2 - aa bR\ 4 t iqi ~ 6t 

We wish to apply this to the symmetrized metric s g a b, which has the prop- 
erty that the rescaled symmetrized metric s g a b satisfies (|6.3[ ). The estimate 
( |6.18| ) translates to an estimate for a metric satisfying ( |6.3| ) after replacing 
ao by ao — 2e, which by the definition of «o satisfies ao — 2e > e > 0. 



Therefore we get in view of the discussion at the beginning of this section 

a , 
d =3 
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or 



Now recall e = ao/4 = min a {p a (xo)}/40 and q\ > 20e by construction. This 
gives 

4 t 3qi 

4 t 3pi 

where we used that fact that q± > pi by construction. 
This finishes the proof of 

Lemma 6.3 (Curvature estimate). 

t 2-a\S R a b ^ t 3pi 

□ 

Now some estimates will be obtained for matter variables. These will be 
used to check that the right hand side of the second reduced system has the 
properties required for a Fuchsian system. Let w a be a one-form with the 
property that w a =4 t 9a ■ We have 

t 2 ~g ab V a w b = t 2 ~g ab ~e a w a - t 2 g ab ~g f9 f abfWg 

Here it has been assumed that w a behaves in a suitable way upon taking 
derivatives. In the context of the matter variables this will be the case for 
the relevant choices of w a , namely e a (cj)) and t~ l u a . Note that this estimate 
requires no use of cancellations, since it only uses the relation ( |6.12a| ) and 
not ( 6.8g| ). In this situation the estimate for a given quantity is never more 



difficult than that for the corresponding velocity dominated part, since the 
difference between the two is always of higher order. This gives the estimates 
required for the matter equations in the scalar field case. For the stiff fluid 
some more work is needed. 

The aim now is to estimate the terms on the right hand side of equations 
( |5.18| ) by a positive power of t. For most of these terms no cancellations 
are required to get the desired estimate. There are only two exceptions to 
this and they will be discussed explicitly now. The first is the following 
combination which arises if the evolution equation for v a is written out 
explicitly: 

(6.19) lT 1 e a (?ii)-ir 1 e a (ji) 
This expression is equal to 

(6.20) t Pl [-A~ l v{\ + A~ 1 ut! 3l )- 1 A- 1 (VaA + V a u) + A'H^Vav) 
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which is 0{t^). The contribution of this expression to the right hand side 
of the evolution equation for v a is as a consequence 0(i^ 1_/32 ) which shows 
that it is necessary to choose 02 < Pi- The second expression where a 
cancellation is necessary is 1 + ttrfc. Now tik = —t~ l + tiKt~ 1+a ° and hence 
1 + tirk = tiKt a ° . It follows that this expression is O(t ao ). 

The analysis of the other terms is rather straightforward, although lengthy, 
and will not be carried out explicitly here. However some comments may 
be useful. The terms which are a priori most difficult to estimate are those 
involving covariant derivatives of u a . For those it is convenient to use the 
components in the rescaled frame e a . For all other terms the original frame 
e a can be used straightforwardly. In order that all terms can be estimated 
by a positive power of t it suffices to choose fi\ and small enough. One 
possible choice is @2 < Pi < Qi — 5e. 

In order to show that the reduced systems for the Einstein-scalar field 
system and the Einstein-stiff fluid systems, are in Fuchsian form, we need 
to show that t 2 ~ aa bM a b = o(t 5 ) for some 5 > 0. 

We consider first the Einstein-scalar field case. In this case, 

Arguing as above for the estimate of t 2 ~ aa b R a b , we have 

Therefore it is enough to show 

t 2 -^M ah ^ & 
By definition M a b = t~ qa ~ Qb e a ((j))eb((j)) and hence 

Mab * 

This give using a a b = \q a - q b \, 

t 2 ~ &a tM a b 4 t 3qi 
For the scalar field case this gives, together with the above, 

t 2-a\ (S R a b _ M aj ^ ^ for gome §>Q 

which is the estimate required for proving that the second reduced system 
for the Einstein-scalar field system is in Fuchsian form. The argument for 
the Einstein-stiff fluid system is very similar. 



7. The constraints 



The main aim of this section is to show that if a solution of the evolution 
equations is given which corresponds to a solution of the velocity dominated 
equations as in Theorem 2.1 or 2.2 then it satisfies the full constraints. It 
will also be shown how the existence of a large class of solutions of the 
velocity dominated constraints can be demonstrated. 
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The first result on the propagation of the constraints relies on rough 
computations which prove the result in the case where all p a are close to 
1/3. An analytic continuation argument then gives the general case. 

Lemma 7.1. Let (°g a b, k ab , (j)) be a solution of the velocity dominated sys- 



tem as in the hypotheses of Theorem 2.1 with \p a — 1/3| < ao/W. Let 



(7%, K a b , ip) be a solution of (5. IS) and ( 5.19Q modelled on this velocity 



dominated solution and define g ab , k ab and (f> by ([3. 1\). Suppose that this 



solution satisfies the properties 1.-6. of the conclusions of the Theorem 2.1. 
Then the Einstein constraints are also satisfied. 

Proof Define: 

(7.1) C = -k ab k ab + (trA;) 2 -R- IQnp 

(7.2) C b = V a k\ - V 6 (trfc) - 8irj b 
These quantities satisfy the evolution equations 

(7.3) d t C - 2(trfe)C = V a C a 

(7.4) d t C a - {tvk)C a = ^V a C 

Define rescaled quantities by C = t 2 ~' ni C and C = t l ~ m C a for some positive 
real numbers rji and 772. Then the above equations can be written in the 
form: 

(7.5) td t C + mC = 2[1 + ttrk]C - t 2 ~ Vl+m V a C a 

(7.6) td t C a + V2 C a = [1 + ttvk]C a - (l/2)^-" 2 V a C 



Choose 771 and 772 so that 771—772 > 0. The aim is to apply Theorem \L2 
to show that C and C a vanish. In order to do this we should show that 
these two quantities vanish as fast as a positive power of t as t — > 0, that 
1+ttrk vanishes like a positive power of t and that the term V a C a is not too 
singular. In obtaining these estimates it is necessary to use the behaviour of 
the derivatives of the solution mentioned in the remark following Theorem 



4.2, Note that since the velocity dominated constraints are satisfied by as- 
sumption, it is enough to estimate the differences of the constraint quantities 
corresponding to the velocity dominated and full solutions, since these are 



in fact equal to C and C a . By property 2. of the conclusions of Theorem 2.1 
it follows that 1 + ttik = O(t ao ) which gives one of the desired statements. 
Similarly it follows that 

(7.7) -k ab k ab + (trk) 2 = -(°k ab )(°k ab ) + (tr°A;) 2 + 0{r 2+a °) 

It follows from the curvature estimates done in section || that the scalar 
curvature is also O(t~ 2+ao ). Now consider the expression p — °p. The 
components of the inverse metric can be estimated by t~ 2+40a ° , so that the 
terms in p involving spatial derivatives can be estimated by t~ 2+a ° as well. 
The difference of the time derivatives can be estimated by t~ 2+ ^ . These are 
the estimates for the Hamiltonian constraint that will be needed. 
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The estimates just carried out were independent of the restriction on 
the p a in the hypotheses of the lemma. The following estimates for the 
momentum constraint are of a cruder type and do use the restriction. First 
note that the metric and its inverse can be estimated by the powers of t equal 
to 2pi and — 2p^ respectively. It follows from the assumption on the p a in the 
hypotheses of the theorem that 2p\ > 2/3 — ao/5 and — 2p^ > — 2/3 — cto/5. 
Using ( |3.4| ) then shows that the connection coefficients can be estimated in 
terms of the power — 2ao/5. The effect on the order of a term of taking a 
divergence can be estimated by the powers — 2ao/5 and —2/3 — 3«o/5 for 
upper and lower indices respectively. 

The gradient of the mean curvature is O(t~ 1+ao logt) while the difference 
of j a is 0(t~ 1+ P log t). The difference of the divergence of the second fun- 
damental form produces the power — 1 + 3«o/5- Evidently the last power 
and that containing (3 are the limiting ones and determine the estimate for 
C a . Similarly the divergence of C a can be estimated by the powers —5/3 
and —5/3 + (3 — 3ao/5. Note that it follows from the definition of ao that 
ao < 1/30. Thus given the hypothesis of the lemma it can be concluded 
that 7/i and r]2 can be chosen in such a way that all terms on the right hand 
side of the propagation equations for the constraint quantities vanish like 
positive powers of t. Thus these equations are Fuchsian and the conclusion 



follows from Theorem 4.2. 



The analogue of this lemma with the scalar field replaced by a stiff fluid is 
also true and can be proved in the same way. Next the restriction on the 
exponents p a will be lifted. Consider a solution (°<? a t>> °&afc; ^>°^a) of the 
velocity dominated constraints for the Einstein-stiff fluid system. Let °k ab 
be the trace-free part of °k ab . The velocity dominated constraints become: 

(7.8) -°k abo k ab + (2/3)(tr°A;) 2 = 16rr/i 

(7.9) V a °k a b = 8tt/iu 6 
Now let x k ab = (1 - Xfkab + (l/3)(tr°k)°g ab and 

(7.10) V = (1/16tt)[-(1 - X) 2 (°k ab )(°k ab ) + (2/3)(tr A;) 2 ] 

(7.11) Vb = 2(1 - X)V a0 k ab [-(l - \) 2 (°k ab )(°k ab ) + (2/3)(tr°A;) 2 ]- 1 

Then °g a b, X k a b, X ^, X Ua is a one parameter family of solutions of the veloc- 
ity dominated constraints which depends analytically on the parameter A. 
There exists a corresponding family of solutions of the velocity dominated 
evolution equations which also depends analytically on A. Next, Theorem 
[4,2| provides a corresponding analytic family of solutions of the full evolu- 
tion equations. (Cf. the second remark following that theorem.) These 
define constraint quantities depending analytically on A. For A close to one 



Lemma 7.1 shows that these quantities are zero. Hence by analyticity they 



are zero for all values of A, including A = 0. This means that the conclusion 



of Lemma 7.1 holds for all positive p a and a stiff fluid. Since any solution of 
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the second reduced system for a scalar field defines a solution of the second 
reduced system for a stiff fluid, this extension also holds for the scalar field. 

A variant of the conformal method for solving the full Einstein constraints 
can be used to analyse the velocity dominated constraints. Consider the 
following set of free data: a Riemannian metric g ab , a symmetric trace- free 
tensor a ab on S and two scalar functions (fi and (fit on S. Next consider the 
following ansatz: 

(7.12) g ab = oj A g ab 

(7.13) k ab = -{l/tyt^gab + uj- 2 l ab 

(7.14) <fi = uj- 2 

(7.15) (fit = uj-% 
where 

(7.16) l ab = a ab + V a W b + V b W a - (2/3)g ab g cd V c W d 



Putting this into ( |2.12a ) and defining p = \{(fit) 2 gives p = uJ~ 8 p, a relation 
well known from the usual conformal method. As a result of the Hamiltonian 
constraint the function (fi satisfies the following algebraic analogue of the 
Lichnerowicz equation: 

(7.17) -^ 12 i ab i cd g ac g bd + |«o 2 " 16™~ 8 P = 

Solving this comes down to looking for positive roots of the equation a£ 3 + 
b^ 2 — c = where a and b are non-negative and c is positive. The derivative 
of the function on the left hand side of this equation is C(^ a C + 26). Thus 
unless a and b are both zero the derivative has no positive roots. Moreover 
the function tends to plus infinity for large £ and is negative at £ = 0. Hence 
the equation has a unique solution for each a and b not both zero and if a 
and b depend analytically on some parameter then the solution does so too. 
If a and b are both zero then of course there is no positive solution. In 
the case of interest here a and b are both positive. The function u) is given 
by uj = £l{l a blcd9 ac 9 bd ,'to:P) where the analytic function O is defined as the 
solution of the algebraic Lichnerowicz equation. The momentum constraint 
implies the elliptic equation 

(7.18) 

g as ^a[V s W b + V b W s - [2/Z)g sb g cd V c W d ] = 8vr[J 6 - 2^V b u] - V a a a b 

for W a . Here j a = (fitS7 a (fi. Note that, when u is expressed in terms of the 
function $7 of the basic variables, it depends on the first derivatives of W a . 
Thus the expression V ' a uj involves second derivatives of W a and is not simply 
a lower order term. 

Consider now the linearization of ( 7.18[ ),with respect to W a , where uj 



has been reexpressed using Q. In particular, consider the linearization in 
the particular case where g ab is the metric of constant negative curvature 
on a compact hyperbolic manifold, the tensor a ab is zero, (fi and (fit are 
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constant and the background value of W a is zero. Because w is a function 
of an expression quadratic in W a , the right hand side of ( [7,18| ) makes no 
contribution to the linearization. Since g a b has no non-trivial conformal 
Killing vectors it follows from the standard theory of the York operator that 
the operator obtained by linearization of the equation (7.18) is invertible 



as a map between appropriate Sobolev spaces. Then an application of the 
implicit function theorem gives solutions of ( 7.18j ) for arbitrary choices of the 



free data sufficiently close (with respect to a Sobolev norm) to the particular 
free data at which the linearization was carried out. This shows the existence 
of solutions of the velocity dominated constraints which are as general as 
the solutions of the full Einstein constraints (at least in the crude sense of 
function counting). 

The conformal method can be applied in a similar way in the stiff fluid 
case and it turns out to be easier than in the scalar field case. This might 
seem paradoxical, since the scalar field problem can be identified with a 
subcase of the stiff fluid problem. The explanation is that it is difficult to 
identify which free data in the procedure for constructing stiff fluid data 
which will be presented correspond to data for a scalar field. The ansatz 
used is \i = lo~ 8 [1 and u a = u 2 u a . This gives the scaling p = lu~ s lo and 
j a = io~ e j a which is often used in the conformal method. The quantities 
describing the geometry are scaled as in the case of the scalar field. The 
equations for lo and W a are very similar in both cases, with the notable 
difference that in the stiff fluid case the term involving the derivative of uj 
is missing from the equation for W a . This means that the equation for W a 
is independent of u and can be solved by standard theory, as long as the 
metric g a b has no conformal Killing vectors. Once this has been done the 
algebraic equation for u can be solved straightforwardly. 

8. Discussion 

We have shown the existence of a family of solutions of the Einstein 
equations coupled to a scalar field or a stiff fluid whose singularity structure 
we can analyse. No symmetry assumptions are made and the solutions are 
general in the sense that they depend on the same number of free functions as 
general initial data for the same system on a regular Cauchy surface. These 
solutions agree with the picture of general spacetime singularities proposed 
by Belinskii, Khalatnikov and Lifshitz in two important ways. Firstly, the 
evolution at different spatial points decouples, in the sense that the solutions 
of the full equations are approximated near the singularity by a solution of a 
system of ordinary differential equations. Secondly there exists a Gaussian 
coordinate system which covers a neighbourhood of the singularity in which 
the singularity is situated at t = 0. It is easily seen that the curvature 
invariant R a ^R al3 = 647r 2 (Va,0V a 0) 2 blows up uniformly for t — > 0. In 
fact the leading term is proportional to A A (x)t~~ A and in the solutions we 
consider A can never vanish, as a consequence of the Hamiltonian constraint. 
Thus these singularities are all consistent with the strong cosmic censorship 
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hypothesis. The mean curvature of the hypersurfaces of constant Gaussian 
time tends uniformly to infinity as t — > so that the the singularity in 
crushing in the sense of JlO[| . It then follows from well-known results that 
a neighbourhood of the singularity can be covered by a foliation consisting 
of constant mean curvature hypersurfaces. This is the most general class of 
spacetimes in which all these suggested properties of general spacetimes have 
been demonstrated. A subclass of these spacetimes is covered by the results 
of Anguige and Tod|l|. The connection between their results and those of 
the present paper deserves to be examined more closely but intuitively their 
spacetimes should correspond to the case where, in our notation, the pt are 
everywhere equal to 1/3. 

The spacetimes constructed have been shown to be general in the sense 
of function counting. It would, however, be desirable to prove that the 
assumption of analyticity of the data can be replaced by smoothness and 
that, this having been done, the spacetimes constructed include all those 
arising from a non-empty open set of initial data on a regular Cauchy surface 
which, in particular, contains the initial data for a Friedmann model. This 
would be a statement on the stability of the Friedmann singularity. A model 
for this kind of generalization is provided by the work of Kichenassamy fl2|l 
on nonlinear wave equations. 

It was indicated in the introduction that the results on the Einstein-scalar 
equations can be interpreted in more than one way. The interpretation which 
has been emphasized here is that where the metric occurring in this system 
is considered to be the physical metric. In the interpretation in terms of 
string cosmology the physical metric is (up to a multiplicative constant) 
e^9iu>- This means that for A(x) sufficiently negative the limit t — > does 
not correspond to a singularity at all, but rather to a phase which lasts for 
an infinite proper time. It is the time reverse of this situation which plays 
a role in the pre-big bang model Q. Another interpretation is in terms of 
the vacuum field equations in Brans-Dicke theory. This is very similar to 
the string cosmology case, with the difference that the conformal factor 
is replaced by e c ^ where C is a constant which depends on the Brans-Dicke 
coupling constant. 

All the results in this paper have concerned the case of three space di- 
mensions. There are reasons to believe that if the space dimension is at 
least ten then the vacuum Einstein equations allow stable quiescent singu- 
larities, similar in some ways to those of the Einstein-scalar field equations 
in three space dimensions |§J. The techniques developed in this paper might 
allow this to be proved rigorously. It would also to be interesting to know 
what happens to the picture when further matter fields are added. There 
are several possibilities here. One is to add some other field, not directly 
coupled to the scalar field, to the Einstein-scalar field system. A second is to 
reinstate some of the extra fields (axion, moduli) which have been discarded 
in passing from the low energy limit of string theory to the Einstein-dilaton 
theory. A third is to add extra matter fields to the Brans-Dicke theory. 
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Another direction in which the results on the Einstein-scalar field and 
Einstein-stiff fluid equations could be generalized is to start with situa- 
tions where the solution has one Kasner exponent negative and investigate 
whether it moves (in the direction towards the singularity) towards the re- 
gion where all Kasner exponents are non-negative. If this were true, then 
the singularities in generic solutions of these equations could be quiescent. 
The set of initial data concerned would be not just open, but also dense. 
This question is sufficiently difficult that it would seem advisable to first try 
and investigate it rigorously in the spatially homogeneous case. 
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